Comment on the recent COMPASS data on the spin structure function gi 
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We examine the recent COMPASS data on the spin structure function g^ singlet. We show that 
it is rather difficult to use the data in the present form in order to draw conclusions on the initial 
parton densities. However, our tentative estimate is that the data better agree with positive rather 
than negative initial gluon densities. 
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I. INTRODUCTION 

The COMPASS collaboration has recently presented new data, Ref. on the singlet component of the spin 
structure function 171 . These data were obtained from the measurements of the longitudinal spin asymmetries in the 
scattering of muons off the LiD target. They found that approximately 

9i{x,Q^)^Q, (1) 

with small errors, in a wide region of x and at small . More precisely, the kinematic region covered in Ref. [H is 

GcoMPAss- 10-" < a; < 10-1; IQ-^ GeV'' < Q'' <l GeV\ (2) 

The fact that 171 is zero in the wide region of x at the first sight looks quite unexpected and even intriguing, and clearly 
requires a theoretical explanation. The most- known theoretical tool for describing 171 is the Standard Approach (SA), 
based on the DGLAP evolution equations, Ref. combined with the standard fits, Refs. [§1 -@, for the initial parton 
densities. However, the values of in the region Gcompass are quite small and therefore this region is beyond the 
reach of DGLAP. In Refs. [1, 0] we have suggested an alternative approach for describing gi at small x and arbitrary 
Q^. Briefly, it combines the total resummation of the leading logarithms of x suggested in Refs. with the shift 

with /i w 5.5 GeV for the singlet gi. This shift automatically leads to an effective change of x: 

X X — X + z (4) 

where z = ^j? /w, w = 2pq, and p, q are the proton and virtual photon momenta respectively. The shifts of and x 
in Eqs. pi4p allowed us to express gi at small x and arbitrary in terms of g^^{x, Q^) obtained in Refs. [9] by the 
total resummation of the leading logarithmic contributions in the region of small x and large Q^: gi at small x and 
arbitrary can be written as gi^{x,Q'^). Let us notice that introducing a shift, similarly to Eq. ([3]), has been a 
common tool for describing the small- kinematic region, see e.g. Refs. and refs. therein. However, contrary to 
all other approaches, we have introduced the shift Eqs. pi4p from the analysis of the Feynman graphs involved, and 
the value of fj, is fixed from theoretical considerations (see Refs. [1,0| for detail). In Ref. p we have predicted that gi 
should not depend on x in the COMPASS kinematic region: indeed Eq. ([4]) shows that gi depends on z rather than 
X in the region Gcompass- The exact value of gi in this region cannot be predicted because it strongly depends on 
the interplay between the quark and gluon contributions. Those contributions involve the coefficient functions and 
the initial quark and gluon densities dq and Sg which are unknown. On the other hand, the very fact that gi = 
approximately, can be used to estimate Sq and 5g. A straightforward application of our results in Refs. [1, 0| to the 
region Gcompass might be misleading. The point is that the approach of Refs. 0] is valid in the kinematic region 
of small X, i.e. for z ^ 1, whereas in the COMPASS experiments 30 GeV^ < w < 270 GeV^ and therefore 



1 ^ ^ ^ 0.1, 



(5) 



2 



In order to extend the approach of Refs. [1, 0| to the region of Eq. ([5]) we suggest in the present paper a simple 
interpolation expression for gi which combines the approach of Refs. [6,, Jj with accounting for the non-logarithmic 
contributions to the coefficient functions in the fixed orders in as- 

The present paper is organized as follows: in Sect. II we remind the expressions for gi obtained in Ref. in 
Sect. Ill we generalize them to the region of Eq. ([5]), adding non-logarithmic contributions to the cocfhcicnt functions 
and anomalous dimensions. Then, in Sect. IV we apply this technique to the COMPASS data. Sect. V is for our 
concluding remarks. 



II. EXPRESSION FOR THE SINGLET gi AT SMALL x AND ARBITRARY IN THE LEADING 

LOGARITHMIC APPROXIMATION 



Explicit expressions for the singlet gi at small x and arbitrary in the Leading Logarithmic Approximation (LLA) 
were obtained in Refs. @, 7|. They account for the total resummation of DL contributions and for the running 
effects. According to Ref. 6|, the LLA expression (= gi^) for the singlet gi at small x and arbitrary is: 
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(6) 



The coefficient functions Cq^\uj) and Cq^' (lu) as well as the exponents D,(^^-^{uj) are expressed through the anomalous 
dimensions Hik which account for the total resummation of DL contributions and the running as effects, as follows. 
We remind here that in our approach we do not use the DGLAP parametrization — as{Q^). Instead, we use the 
alternative we suggested in Ref. [l^l . It allows us to consider the region of really small Q'^ and at the same time to 
be within the framework of the Perturbative QCD. 



A. Expressions for the exponents fi(±) 
The explicit expressions for f^(±) are: 

n(^±)^l[Hqq+Hgg±VR\ , 



where 



R - [Hqq - Hggf + AHqgHgq 



and Hqq, Hqq, Hqq, Hqq are the anomalous dimensions calculated in LLA. 



B. Coefficient functions 



The expressions for the coefficient functions are also written in terms of Hn 
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Here 



^ ^ ^99 ^ ^qq- T — Uj'^ ~ ^{Hgg + Hqq) + {HggHqq — HgqHqg) 



(10) 
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C. Anomalous dimensions 



Here we have: 



H. 
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(11) 



where 



Z 



V2 



bik = aik + Vik 



(12) 
(13) 



with the Born contributions aik defined as foUows 
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with 77 = ln{p^/A^Q(j^) and 6 = (33 - 2n/)/(127r). 



p + rj 



1 
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(19) 



III. EXPRESSION FOR THE SINGLET gi AT ARBITRARY x AND 

Our goal now is to obtain exphcit expressions for the singlet gi which could be valid at arbitrary x and . The point 
is that the Eqs. ([9|) and (jlip for the coefficient functions and anomalous dimensions present the total resummation 
of the leading logarithms of x but those contributions are large when x <C 1 only. Alternatively, non-logarithmic 
contributions can be large at large x 1 and should be taken into account at large x. Such terms are beyond the 
rich of our approach, so we cannot do the total resummation of them. Instead, we can obtain them in the orders 
^ Us and ~ a^. Adding these contributions to the expressions in Eqs. ([9|) and ()11|) . we arrive at new formulae for the 
coefficient functions and anomalous dimensions, which are valid at arbitrary x. In doing so, we can use the DGLAP 
results for the anomalous dimensions and coefhcient functions. Let us demonstrate it in detail, using an example of 
the singlet coefficient function C,. The dealing with the other coefficient function and anomalous dimensions is quite 
similar. The NLO DGLAP singlet coefficient function C^'^^^^ in the lu -space and at integer lu = n is (see e.g. |11| ) 



C, 



DGLAP 



1 



27r 



1 



7i{n + 1) 



2n 



1 



(20) 
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where we use the standard notations 

n n 

Si in) = ^ 1/fc , ^2 (n) - 5] 1/fc^ (21) 

k=l k=l 

The expression ([20)1 is obtained by direct calculation of the Feynman graphs and is insensitive to the value of Q^, save 
the parametrization of a^. So, we can borrow it, though after some appropriate changes: In the first place it should 
be valid at arbitrary w; second, according to the results of Ref. 13], the coupling as{Q^) should be changed to A{Ld) 
defined in Eq. (fT7|) . The analytic continuation of Eq. ([20|) to arbitrary uj is obtained through expressing the sums in 
Eq. (|2ip in terms of the polygamma ■0 -function and the Euler constant C: 

Si{n) = C + n^{n-l), S2{n - I) = — + i,' {n) . (22) 
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After that we obtain an expression which we address as clp accounting for both logarithmic and non-logarithmic 
contributions in the first loop. Repeating the same procedure for the gluon coefficient function, we obtain its first-loop 
value clp. Apart from the trivial replacement 81,82 by '4^{lo) according to Eq. (P^ . Cg^-* and C^^^ differ from the 
NLO DGLAP coefficient functions C^^° ^glapj-^) ^^ve treatment of a,: 

with A being defined in Eq. p7|) . The two- loop expressions h\^^ for the anomalous dimensions can be found quite 
similarly. They also can be obtained from the NLO DGLAP anomalous dimensions ^'-^^^^{oj) with expressing 

Si, 82 through ^piuj) and replacing as{Q^) by A{uj): 

i?L'H-) = 7ir"°^'^"^^(-)U.-.A. (24) 

Explicit expressions for the NLO DGLAP coefficient functions and anomalous dimensions can be found e.g. in Ref. [ll|. 
Obviously, the replacement as{Q^) by A{uj) in Eqs. (|23l24p makes possible to use C^^\C^^^ and hI^^ at arbitrary 
in contrast to the DGLAP expressions for the coefficient functions and anomalous dimensions. Combining Cq^"* , Cg^^ 
and H^^^ with Eqs. (|9llip . we obtain the interpolation formulae equally valid for small and large x. Indeed, the 
replacements Hik by Hik and Cq^^ by Cq^^ in Eq. ([6]) allow to extend the small- 5; formula Eq. ([6]) to arbitrary x. 
The new coefficient functions Cq^^ are defined as follows (the superscripts ± are dropped here): 

Cq^Cq+ - ACq , Cg=Cg + C^^^ - ACg (25) 

where Cq,g are defined in Eq. ACq^g are their perturbative first-loop expansions and C^^\ Cf^ are given by 
Eq. (I23p . The definitions for new anomalous dimensions Hn^ look quite similar: 

H,k = H,k + h\1^ - AH,k (26) 

where Hik are introduced in Eq. ([TT|) . AHik include the first and second terms of their perturbative expansions 
whereas H^^, are given by Eq. (|24p . Now, introducing f^(±) according to Eq. ([7]), with Hik in place of Hik, we arrive 
at the expression describing gi at arbitrary x and Q^: 

o2 
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When fi^, Eq. ((27l) can be expanded in the series in Q^/fi^ 
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9i{x,Q^) « gi{z) + (gV;.2)Mjl^ + 0{{QV^,'f) (28) 
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where 



We have denoted here 



and 
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(29) 



(30) 



, AC, 



(31) 



IV. IMPLICATIONS FOR THE RECENT COMPASS DATA. 



Eq. ((28|) shows exphcitly that gi practically does not depend on x in the region of Eq. ([2]). It perfectly agrees with 
the flat a;-dependence of gi observed experimentally in Ref. Such a dependence means that gi in the COMPASS 
kinematic region ^ does not depend on the conventional variables x and . On the contrary, Eq. predicts 
that the z-dependence of gi is pretty far from being trivial. Let us notice here that z is inversely proportional to the 
standard variable v — w/{2M) measured in GeV, with M — 1 GeV: 



\2M/ V 



so the region ([5]) covered in the COMPASS experiment corresponds to the z/-region (in GeV) 



15 < 1/ < 150. 



(32) 



(33) 



Obviously, a straightforward and unambiguous application of our description of gi to the COMPASS experiment could 
be obtained just by fitting the COMPASS data on gi{z). Unfortunately, this is impossible because the COMPASS 
collaboration has not studied the z-dependence of 171. Nevertheless, it is clear that Eq. ([1]) could be satisfied at any z in 
the region ^ only if there exists a strong correlation between 5q and 5g to compensate the difference between Cg and 
Cg explicitly given in Eq. ([?!]) . We think that the chance for such a correlation is very tiny, though strictly speaking 
this situation cannot be excluded. An alternative interpretation of the COMPASS result is to consider Eq. ([T]) as 



< gi{z) >= 0. 



(34) 



where < gi{z) > is the average value of gi observed by COMPASS. Obviously, in order to match Eq. gi{z) should 
acquire both positive and negative values in the region ([5]) . For further investigations with Eq. (I34p one should choose 
appropriate fits for the initial parton densities Sq{z) and Sg(z). Such fits are practically absent in the literature. In 
Ref. [6| we suggested to approximate Sq{z) and 5g{z) at small z by constants to get a rough estimate. However, z in 
the COMPASS region ^ is not small, so we prefer to use a DGLAP-like set of fits: 



Sq{z) = Ngz{l - zf{l + 3z), dg{z) = Ng{l - z)'*(l + 3z). 



(35) 



This set corresponds to the DGLAP-fits suggested in Ref. ^ but does not coincide with them. The difference is in 
the power factors z°- while the terms in the brackets in Eq. (|35p and in Ref. Q coincide ( a; in Ref. @ is replaced by z 
in Eq. ((35)) ). Indeed the fit for 5q in Ref. [3] contains the singular power factor x~^^^ whereas the power factor for Sg 
is x'^-^. In Ref. 12] we have proved that the role played by the singular terms x~°' in the DGLAP fits is to mimic the 
total resummation of ln'^(l/a;) . When the resummation is taken into account, such factors ( namely the factor z~'^'^ 
in Sg) does not make sense any longer and should be dropped. The same is obviously true when x is changed by z. 
So, extracting the singular factor a;^" from the fit in Ref. Q, we arrive at Eq. ([55]) . Now it is easy to check that the 
fits ([55]) do not lead to a flat z-dependence for gi and cannot keep gi{z) = in the whole COMPASS region 

In more detail by substitution of Eq. (j35p into Eq. (|29p and performing the integration over uj numerically, with 
fixed and positive Nq and varying the values of Ng, we plot our results in Fig. 1. By a close inspection of the various 
configurations shown, we can easily conclude that these fits could be compatible with Eq. only if A^g > and 

Ng > Ng. 
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As the way of averaging gi over z in the COMPASS data is unknown, we can try another possibihty, approximating 

<gi(z)>«gi(<z>) = 0, (36) 

where < z >= 0.25 (i.e. < i/ >« 60 GeV) is the mean value of z from the region ((5]). Then using Eqs. (j29|35|) . 
keeping positive Ng and varying Ng. Figs. 1 suggest again that Ng are positive and Ng > Ng. 

V. CONCLUSION 

In the present paper we have considered in detail the recent COMPASS data on gi. These data first confirm our 
prediction in Ref. ^ that gi at small does not depend on and x. Instead, we predict that gi depends on the 
invariant energy w — 2pq and the experimental investigation of this dependence would allow to estimate the initial 
parton densities. Unfortunately, this information is absent in the present COMPASS data, so a reliable study of the 
initial parton densities cannot be done. However, we have suggested two possible interpretations, Eqs. (|34p and (|36p . 
of the COMPASS result Eq. ([1]). Combining the LLA resummation with the explicit first-loop values of the coefficient 
functions and using the DGLAP-like parametrization ([55]) of the initial parton densities, we conclude that the data 
suggest rather positive than negative values of the initial gluon density. We remind that our analysis is tentative. 
More quantitative conclusions can be drawn only after an accurate experimental study of the z-dependence of gi has 
been performed. 
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FIG. 1: The v -dependence of g\(y\ with &q, 5g defined in Eq. (|35|l . for Nq = 0.5 and different values of Ng-. (a) -1.5, (b) -0.5, 
(c) 0, (d) 0.5, (e) 2, (g) 3.5; the COMPASS v -region corresponds to Eq. ^ . 



